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COHERENT  AND  INCOHERENT  RADIATION  FROM  FREE  ELECTRON  LASERS 

WITH  AN  AXIAL  GUIDE  FIELD 

I .  INTRODUCTION 

Recent  experiments^""^  have  amply  demonstrated  the  usefulness  of 

the  free  electron  laser  as  a  tunable  source  of  coherent  submillimeter 

radiation,  and  have  stimulated  interest  in  a  variety  of  experimental 

configurations.  Theoretical  work^  °  in  this  regard  has  concentrated  on 

the  case  in  which  a  relatively  low  density  (i.e.,  self-field  effects  are 

negligible)  relativistic  electron  beam  is  propagated  through  a  periodic 

helically  symnetric  magnetic  field,  referred  to  as  the  "wiggler"  field. 

In  the  beam  frame  the  wiggler  field  appears  as  a  backwards  propagating 

electromagnetic  wave,  and  radiation  is  produced  by  means  of  coherent 

scattering  of  this  wave  off  electromagnetic  fluctuations  in  the  beam. 

Another  configuration  employed  is  one  in  wliich  an  axial  guide  field  is 

present,  and  is  used  principally,  but  not  exclusively,  in  the  regime  in 

which  electron  densities  and  currents  are  sufficiently  high  that  the 

axial  field  is  required  to  contain  the  beam  against  the  effects  of  the 

/ 

self-fields.  Theoretical  analyses  of  the  coherent  radiation  mechanism 

in  this  regime  have  begun  to  appear  in  the  literature  for  this  config- 
,,  7-11 

uration  as  well. 

It  is  our  purpose  in  this  work  to  treat  both  the  coherent  and 
incoherent  radiation  from  a  free  electron  laser  configuration  which 
contains  an  axial  guide  field,  in  the  low  density  limit.  Full  relativ¬ 
istic  effects  are  included,  but  the  self- fields  of  the  beam  are  neglected. 

The  principal  application  of  this  study  to  existing  experiments  is  the 

1  n 

Stanford  University  free  electron  laser  *"  in  which  a  low  current  (~  2  A), 
high  energy  (24-43  MeV)  beam  was  injected  into  a  drift  tube  in  which 
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an  axial  guide  field  of  1  kG  and  a  wiggler  field  of  about  2.4  kG  were 
applied.  In  this  regime  collective  beam  effects  are  negligible. 


The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we 
discuss  the  single -particle  trajectories  of  electrons  in  a  magnetic 
field  geometry  conposed  of  a  combined  axial  guide  field  and  helical 
wiggler  field.  The  incoherent,  spontaneous  emission  spectrum  is  cal¬ 
culated  in  Sec.  Ill  using  test  particle  techniques,  and  the  linear  gain 
is  found  in  Sec.  IV  in  the  limit  of  a  low  density  beam.  In  treating  a 
low  density  beam,  we  restrict  the  analysis  of  the  linear  gain  to  the 
small- signal  regime  but  provide  a  fully  kinetic  derivation  of  the  gain 
based  upon  the  Vlasov-Maxwell  equations.  This  is  in  contrast  to  the  work 

n  I  A 

of  Kwan  and  Dawson  and  Bernstein  and  Friedland  in  which  collective  eff¬ 
ects  due  to  high  beam  density  were  included  in  the  context  of  a  fluid  theory 
of  the  interaction.  The  small-signal  gain  has  also  been  considered  by 
Friedland  and  Hirschfield^'''  by  means  of  a  fluid  analysis.  Thus,  we  provide 
a  fully  kinetic  expression  for  the  small-signal  gain  and  an  extensive  des¬ 
cription  of  the  parametric  dependence  of  the  gain  in  the  cold  beam  limit. 

A  summary  and  discussion  is  given  in  Sec.  V. 
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II. 


SINGLE-PARTICLE  TRAJECTORIES 


The  wiggler  field  is  generally  due  to  a  helical  current  winding 
and  the  resultant  field  can  be  shown,  in  a  self-consistent  way,  to  be  of 
the  fomi12,lJ 

&'  *Vr’z)>  (1) 

where  the  wiggler  field  is  derived  from  a  vector  potential  of  the  form 

=  *  TT  (ft  X1  ®wr)  C0S  (e  '  kwZ)  ~r 

w  \  w 

-  I.  (k  r)  sin  (9  -  k  z)  e.  J  m 

1  w  w  ~9 J  ,  [Zj 

in  cylindrical  coordinates.  In  Eq.  (2),  Ew  and  k*  (=  where  ^  is 

the  wiggler  period)  are  assumed  to  be  constant,  and  1^  and  1^  are  the 
modified  Bessel  function  of  the  first  kind  and  its  derivative  respectively. 
In  most  free  electron  laser  experiments,  however,  the  initial  beam  radius 
is  a  small  fraction  of  the  wiggler  period,  and  expansion  in  powers  of  l^r 
yields 

A  A 

B  =-  B,,  (e  cos  k  z  +  e  ,  sin  k,2) .  (3) 

~w  w  w  ~y  * 

Rather  than  work  in  the  laboratory  frame,  we  choose  to  transform 
to  the  frame  rotating  with  the  wiggler  field.  To  this  end  we  define 

A  A  a  aa  a  aa 

e.  =  e  cos  k  2  +  e  sin  k  z,  e,  ■  -  e„  sin  k,  z  +  e,r  cos  k  z,  ex  =  e_  and 

~1  ~x  w  ~y  w  ~Z  ~x  w  ~y  w  ’  ~3  ~z 

write  the  orbit  equations  corresponding  to  the  field  structure  given  in 

(1)  and  (3)  in  the  form 
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(4) 


41  =  v2  (kwv3  •  V  ’ 

*2  ‘  '  Vs  '  Vl(kwv3  '  V’ 

V,  =  Q  V7  , 
i  W  2  ’ 

2  2 

where  ^  =  |eBg/ymc|,  =  leB^ymcl  ,  and  y  =  (1  -  v  /c  )  .  Since  the 

2  _  2  2  2 

total  energy  is  a  conserved  quantity  v  =  v^“  +  v^  +  is  a  constant 
of  the  motion  and  Eqs.  (4)  are  fully  relativistic.  The  principal  benefit 
derived  from  this  transformation  is  that  the  equations  of  motion  are  now 
coordinate -independent,  and  depend  only  on  the  components  of  the  velocity 
and  their  derivatives. 


In  addition  to  the  total  energy,  a  second  constant  of  the  motion 

can  be  found  by  elimination  of  v?  from  (4);  specifically,  u  =  v1  ■  kw  x 

(v3  "  ^  /kw-^2V  is  analogous  to  the  axial  invariant  discussed  by 

14 

Davidson  and  Uhm,  which  in  the  B  .  ■*  0  limit  implies  that  the  axial 

w 

velocity  is  conserved  and  defines  the  rotating  frame  to  be  that  of  the 
Larmor  rotation  of  an  electron  in  a  uniform  magnetic  field.  By  means  of 
these  constants,  the  problem  can  be  reduced  to  the  solution  of  a  single 
nonlinear  differential  equation  for  v, , 

(dx/dt)2  +  <S>  (x)  =  0,  (5) 

where  x  =  6-  -  $  x  =  ft  t/2,  6,  e  v./c,  8  =  !)  /kc, 

u  00  0  Q  w 

*(X)  =  x4  +  4£80(£80  +  8U)  x2  +  8e2BQ3x  +  4e2802(302  +  8U2  -  8/), 


e  =  B^/Bg,  3y  =  v/c,  and  8u  2  u/c.  Note  that  Eq.  (5)  has  the  trivial 
solution  x  =  0  (i.e.,  v_  =  constant)  in  the  limit  of  a  uniform  axial 
magnetic  field  (B  -*•  0) . 

A  detailed  discussion  of  the  solutions  of  Eq.  (5)  has  been  given 
by  Freund  and  Drobot;'*'3  however,  we  restrict  our  attention  here  to  solutions 


4 


corresponding  to  relatively  uniform  axial  velocities.  The  reason  for  this 
is  that  the  radiation  mechanism  is  a  resonant  one  in  which  the  emission 
frequency  is  given  by  u>  *■  2yz  kwVj,  where  yz  =  (1  -  v^  /c  )  2 .  As  a  result, 
variations  in  the  axial  velocity  of  the  order  of  Av^  lead  to  a  broadening 
of  the  emission  spectrum  which  scales  as  Acd  ~  2y7\wAv3.  Thus,  while  small 
oscillations  in  v^  about  some  bulk  axial  velocity  lead  to  a  relatively 
narrow  bandwidth,  large  oscillations  can  result  in  spectral  broadening  with 
a  corresponding  decrease  in  the  linear  gain. 

The  conditions  which  lead  to  small  oscillations  in  the  axial 
velocity  can  be  investigated  by  consideration  of  the  roots  of  the 
pseudopotential  <s(x) .  Physically  meaningful  (i.e.,  real)  solutions  are 
possible  only  when  $  <  0;  therefore,  the  real  roots  of  the  pseudopotential 
correspond  to  the  bounds  on  the  oscillation  about  some  bulk  axial  velo¬ 
city  (defined  by  the  local  minimum  in  4>  between  the  roots) .  The  typical 
character  of  the  pseudopotential  when  e  <  1  and  Sg  >  1  is  shown  in  Fig. 

1.  In  this  regime,  the  pseudopotential  has  two  real  roots  which  for 
given  values  of  e ,  3q ,  and  v  shift  in  both  position  and  spacing  with 
variations  in  6u  within  some  fixed  range.  The  extrema  in  the  range  of 
3u  occur  when  the  real  roots  of  $  are  degenerate,  and  correspond  to 
solutions  with  constant  axial  velocity^ 
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2  2? 

where  the  axial  velocity  v;j  is  given  by  v  =  v1  +  v 3“  for  some  choice 
of  the  total  energy. 

The  values  of  C=  BQ  +  x)  corresponding  to  the  real  roots  of 
$  are  shown  in  Fig.  2  versus  Su  for  e  =  .1  and  B^  =  1,  and  4.  Thus,  the 
figure  describes  the  bounds  on  the  oscillations  in  the  axial  velocity  for 
given  Bu-  It  is  important  to  recognize  that  small  deviations  in  8U  from 
the  extrema  result  in  trajectories  which  do  not  differ  greatly  from  those 
with  uniform  axial  velocities.  Because  of  this  two  electrons  which  are 
characterized  by  only  slightly  different  values  of  8u  and  which  are  close 
together  at  some  point  will  not  undergo  a  large  separation  in  the  course 
of  their  orbits.  To  be  more  specific,  an  initially  bunched  electron  beam 
will  tend  to  remain  bunched.  This  describes  an  "orbitallv  stable" 
uniform-v.  solution,  in  the  sense  discussed  by  Friedland.^  A  quantitative 
test  for  orbital  stability  can  be  expressed  in  the  form 


which  (since  Vm  <  c)  is  trivially  satisfied  for  the  uniform-v-  trajectories 
when  Bq  >  1. 

In  contrast,  when  Bq  is  less  than  unity  the  pseudopotential  can 

have  as  many  as  four  real  roots  and  as  many  as  four  uniform-v.  orbits 

appear.  An  example  of  this  is  shown  in  Fig.  5,  where  we  plot  6-  corresponding 

to  the  roots  versus  3  for  e  =  .1,  and  {L  =  *5  and  1.  The  orbits  with 

u  0 

uniform  axial  velocity  correspond  to  the  points  of  vertical  slope  in  the 
figure,  and  when  Sq  =  .5  four  such  trajectories  appear.  Three  of  these  orbits 
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are  stable;  however,  the  trajectory  corresponding  to  point  A  in  the  figure 
represents  a  fundamentally  different  class  of  solution  than  the  other 
three.  The  behavior  of  the  pseudopotential  for  values  of  Bu  in  the 
vicinity  of  a  stable  orbit  is  shown  in  Fig.  1.  In  contrast,  the  behavior 
of  $  with  Bu  in  the  vicinity  of  point  A  is  shown  in  Fig.  4  in  which  the 
uniform-v-  trajectory  is  obtained  when  the  central  maximum  is  zero.  Arbi¬ 
trarily  small  deviations  from  this  condition  result  in  drastically  different 
types  of  trajectory  with  large  fluctuations  in  the  axial  velocity.  Such 
uniform-v-  trajectories  are  referred  to  as  being  "orbitally  unstable," 
since  a  bunched  electron  beam  with  a  spread  in  the  axial  invariant  in  the 
vicinity  of  point  A  will  rapidly  disperse. 

.As  Bq  decreases  further,  the  range  of  3u  in  which  $  possesses  four 

real  roots  increases,  as  shown  in  Fig.  5  for  Bq  =  .1.  In  particular,  the 

high  axial  velocity,  uniform-v-  trajectory  for  motion  parallel  to  Bq  becomes 

more  pronounced  as  Bq  decreases  further  relative  to  e.  This  characteristic 

of  the  solutions  in  the  small-Sp  regime  exists  even  where  e  >  1  (see  Fig.  6). 

Finally,  we  conclude  that  the  requirement  that  electrons  propagate  with 

relatively  small  fluctuations  in  axial  velocity  is  most  difficult  to  satisfy 

when  Bq  ^  1.  The  natural  corollary,  therefore ,  that  extremely  broadbanded 

17  18 

emission  will  occur  in  this  regime,  is  supported  by  observation.  ’ 

In  order  to  treat  cases  of  nearly  uniform  axial  velocity  we 
implicitly  restrict  consideration  to  trajectories  corresponding  to  the 
stable  uniform-v^  solutions.  In  this  limit,  we  can  write  approximate 


solutions  for  the  momenta  in  the  laboratory  frame  as 


px  2„  -  kv 
0  w  ii 


cos  k  z  +  P  cos  n_t  -  P  sin  fi-t, 
w  x  0  y  O’ 


py  =  nn  -'k  v  sin  V  +  Px  sin  Q0Z  +  Py  cos  *0*' 

U  w  ii  J 


pz =  pn  -  gr-.-is  v  L?x  cos  (V  •  n0t}  •  py  sin  ■  not)] 

u  w  ■*  ' 


where  p(|  ,  P^  and  P^  are  approximate  constants  of  the  motion,  and  vM  = 

pM/Ym.  Observe  that  in  the  limit  of  a  vanishing  axial  guide  field  Eqs. 

(7)  recover  the  well-known  result  in  which  P  and  P  are  the  canonical 

x  y 

momenta  in  the  transverse  direction.  This  approximation  remains  valid 
2  ?  2 

as  long  as  P^  +  P^“  «  p(|  ,  and  the  total  momentum  is  given  by 


to  within  terms  of  order  c2.  The  form  of  these  trajectories  is  shown 
schematically  in  Fig.  7.  Particle  motion  predominantly  follows  the 
helical  wiggler,  but  also  contains  a  small  component  of  Larmor  procession 
due  to  the  axial  field. 
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III. 


SPONTANEOUS  EMISSION 


The  time-averaged  radiated  power  P  can  be  computed  by  means  of 
the  equation 


T/2 

P  =  -  lim  i-  f  dt  f  d3x  JgCfct)  •  , 

T~°  -1/2 


(9) 


where  ^(jc,t)  denotes  the  microscopic  radiation  field,  and 

Nb 

JC^c.b)  =  -e£^jW  6(!C-  SjOO)  (10) 

j=l 

is  the  source  current  composed  of  the  sum  of  the  microscopic  currents 
due  to  all  N^  electrons  in  the  beam  during  the  interaction  time 
-  T/2  <t  <  T/2.  The  radiated  power  can  be  expressed  in  terms  of  the 
Fourier  amplitudes  of  the  microscopic  fields  and  source  currents  in  the 
following  manner 

00 

P  =  -  2(2ir)4  lim  £  / d3k  / dw  Re  (E,  •  J *  ) ,  (11) 

T-*00  o  v  V 

where  the  asterisk  (*)  denotes  the  complex  conjugate,  and  the  Fourier 

-f  00 

amplitude  is  defined  as  f,  =  /  d4  /  du>  exp  (icut  -  ik  •  x)  f(x,t). 

A  self-consistent  relation  between  the  fields  and  source  currents  de¬ 
pends  upon  the  dielectric  properties  of  the  beam,  and  can  be  written  as 


% 


4Tri  T 

-  «  4, 


A,  • 


(12) 


where  the  dispersion  tensor 


7 


it,*'  (1  •n') 


(13) 


n  (  =  ck/  is  the  index  of  refraction,  £  is  the  unit  dyadic, 


% 
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£  e  £/|£| ,  and  e,  is  the  beam  dielectric  tensor.  For  the  case  of  a 
diffuse  beam,  in  which  the  radiation  frequency  greatly  exceeds  the  electron 


plasma  frequency,  c,  -  ^1  and  (12)  can  be  inverted  to  give 


-  (4rrij., )  / 


*,«,  •  n2W  '  &  jJ  • 

L  %  'A, 


(14) 


As  a  result,  if  the  emissivity  n(o),  ft£)  is  defined  to  be  the  power  radiated 

per  unit  frequency  and  volume  into  the  solid  angle  subtended  by  k, 

a. 

then  we  find 

2 


nfc»,  AT)  =  (27T)°  (*2/Vc3)  lim  1  (|J 


"k 


T-h» 


k,u) 


^  *  ^k,  w 
% 


2 

I  ) 


k  =  w/c 


,  (15) 


where  V  is  the  total  volume  of  the  interaction  region.  The  radiation 

spectrum,  therefore,  can  be  determined  from  (15)  with  a  knowledge  of  the 

single-particle  orbits  (7)  necessary  to  compute  the  source  current. 

Using  the  trajectories  in  (7),  we  find  that  after  transforming 

to  the  frame  (e,  =  cos  \p  e  +sin  ^  e  e,  =  -  sin  ^  e  +  cos  ip  e  ,  and 

v  'v y  ^X  r^y 

.  fk  =  /  k7  +  k2  1  i 

Z  r\j  Z 


tan  -p  =  k /k  )  in  which  k  =  k,  e  ,+  k  e  fk  =  / k 7  +  k2  )  the  source 

y  X  Oj  i  n.  A  Z  o.Z  X  V 


y 


current  is 


Jv 


lim  ^ 


3  t—t  i  - 133;;;'''  J)  ^  imn 

(2tt) ^  T-oo  j=i  0  £ ,m,n= 


sin [(u>-k,vP ^ - (£+n)k  vP ^ - (m-n)Jl^ ^ )T/2]  m  m  m 

"  11  w  11  -2 - -  exp[ib1UJsin(i>UJ-4))-im((})Uj-4>)J 


tt (oj-k,\|p  *  -  (£+n) lTy,  Ll  -  -  (m-n) ^ ) 


(16) 


x  exp[ib^ ^ sin(kw:^  ^ -ij>) -iJtfk^-1 ) -i|»)] exp [i^J ^sin(ku,tp ^  -<p ^ ^ )  -in(kw=p  ^  -$ ^ )] 
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- : 


where  defines  the  initial  position  of  the  j^1  electron, 

<J>^  =  tan"1(Py/PJC)t  b^  =  kx^Vkw(n^-kwV|j^) ,  b[^  =  k^pVfipP 
b^  =  k7Vp-*^ V(S1^ -kwv||^) ^,and  =  (P^  +  P“)/ym.  In  addition, 


y(j)  =  j  (b^Xl 
-Uim  nl  z  J 


-) 


jtk  i  r -s 

"  1  JI(b(^)J„(l>P))^ 


*  ‘(“-w  )jI(b°))Jm(bp))  *  Vp’j^bOb^Cbi'1’))^ 

♦  (vp)  .  nJHt 


>(j).V  v(j) 

w 


ki 


’k"v" 


(17) 


where  vw  =  I^Vjj  /(Q^-k^  ),  and  and  are  the  regular  Bessel  function 
of  the  first  kind  and  its  derivative  respectively.  In  computing  the 
quadratic  forms  of  the  source  current  which  appear  in  the  emissivity  (15) , 
we  impose  a  random  phase  approximation  to  obtain 


n(w, 


¥ 


-  e  CJ 


V  oo  CO 

IE  E 


Zttc^V  j=l  2.,m,n=  -oo  i'm' ,n'=- 


^JZ+n.JZZ+n  ^m-n,m'-n' 


x  l(vj 

V'vZmn 


-^.VpP!  k-Vp^-vJj^Ww-k  vf^-(l+n)k  ^^-(m-n)fl^- 

ran  'bZ  ra  n  K  'v-Jlmn^  m  rr  zj|  ^  J  wj|  v  q 


(18) 


k=^/c 


We  treat  the  emissivity  in  the  limit  in  which  bp^  <<  1,  which  is  gen- 
erally  valid  as  long  as  2y  ■£Bq(V1/v||  )  <<  1.  When  this  condition  is 
satisfied  the  dominant  terms  in  (18)  are  those  for  which  n  =  n'  -  0, 
and  the  emissivity  takes  on  the  relatively  simple  form 


11 


where  0=  cos  ^ (k-£.J  is  the  polar  angle  between  the  wave  vector  and  the 
axial  guide  field. 

We  now  convert  the  discrete  sum  over  individual  electrons  in  (18) 
and  (19)  into  a  continuous  integral  over  the  beam  distribution  function 


by  making  the  replacement  V"  1  5*  ^  / dpxdpydPFb  >  where  nb  =  -\/v  and 
tlie  beam  distribution  is  assumed  to  be  a  function  of  (Px,Pv,p).  As 
a  consequence, 

2 


O  1 

rn  ^  r 

>  — V  JdPxdPvdPFb 


T-l 


7TC 


Yh  cos0-''j|  sine] 


l 


cos6-\|j  sinej  J“  Cb)  J^(b 


x  6(u>-k,v(,  -Ik  v, 


k=o)/ c 


(20) 


where  v..  is  determined  by  means  of  Eq.  (8).  Since  (8)  is,  in  general, 
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II 


a  quartic  polynomial  for  vj|  ,  care  must  be  exercised  in  selecting  the 
appropriate  root.  In  the  limit  of  propagation  parallel  to  the  axial  guide 
field,  (2.0)  reduces  to  the  comparatively  simple  form 


2  2 
e  n^ui 


;  j: 

AS  'vz  2ttc 


.  (v  2 

JdPxdPydpFb  j-V  6  (a)-  (k+kw)  Vj| ) 


+  — T  5(cj~ky,  -ft  )) 

c”  II  0  )  k-u/c 


(21) 


Thus,  the  incoherent  radiation  corresponding  to  the  usual  free  electron 

?  2 

laser  resonance  at  oj  ^  2y“k  v M  scales  as  v“.  The  emission  at  the 

I  W  ||  W 

cyclotron  resonance  is  expected  to  be  much  less  intense  in  this  regime 

■> 

in  which  it  is  required  that  Vt  <  v“  in  the  derivation  of  the  relativistic 

trajectories  (7).  .An  important  consequence  of  this  analysis  is  that  the 

emissivity  in  the  presence  of  an  axial  guide  field  is  enhanced  by  a 

_2 

factor  of  (ft„/k  v  -  1)  .  Note,  however,  that  there  is  no  singularity 

U  w 

2  2  2  2 
since  v  must  be  computed  from  vi  +  v  =  v  <  c  . 
li  it  w 

In  order  to  illustrate  typical  noise  levels  and  spectra  to  be 
expected  in  the  course  of  operation  of  a  free  electron  laser  we  assume  a 
distribution  of  the  form  F(P  ,Pv,p)  =  nb£fPx)ifPy)CbCp) ,  and  model  the 
variation  in  total  momentum  with  a  distribution  of  the  form 


VP) 


= 

17 


tp-Po)' 


ip 


where  p^  characterizes  the  bulk  momentum  in  the  beam  and  ip  describes  the 
momentum  spread.  Note  that  since  P  =  Pv  =  0  the  component  of  the  motion 
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describing  Larmor  rotation  is  assumed  to  be  negligible,  and  the  noise 
spectrum  is  predominantly  associated  with  the  free  electron  laser  resonance. 
As  a  consequence,  the  emissivity  is  of  the  form 


nC“,°fc4)  ~  7373 


U  .2  ml  4  (lnA  )3/“ 

^  W  _ b  u) _ pn _ 


Au 


^  (-+kwc^3  ^2o^o *V2  (uPh'uo)Z  + 


2 - 2  »  (22) 


Au 


?  2 

where  =  47re  n^/m,  Au  =  Ap/mc,  Ug  =  pg/mc,  and 


2  „  2  2 
1  *  w 


ph  Fc  Y  !  +  2w/kwc  • 


As  an  example,  we  shall  evaluate  the  emissivity  (22)  for  a  1  Amp 

electron  beam  characterized  by  a  bulk  energy  of  1.4  MeV  with  a  3%  energy 

spread  and  a  3  cm  beam  radius.  The  beam  density,  therefore,  is  approximately 
9  -3 

n^  i.  1  x  10  cm  .  In  addition,  we  shall  assume  the  wiggler  amplitude  and 
period  to  be  700  G  and  3  cm  respectively,  and  the  axial  guide  field  to  be 
15  KG.  The  resulting  noise  spectrum  is  shown  in  Fig.  8.  The  peak  intensity 
is  of  the  order  of  .013  pW-Hz  ^-cm  3  and  occurs  at  a  frequency  of  approximately 
240  GKz  with  a  half -width  of  6?*. 
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IV. 


COHERENT  EMISSION 


In  this  section  we  derive  the  linearized  gain  by  solution  of 
the  Vlasov- Maxwell  equations.  If  the  distribution  is  written  as  the 
sum  of  equilibrium  and  fluctuating  components  + 

<-f^(z,^,t) ,  then  the  formal  solution  of  the  Vlasov  equation  for  the 
perturbed  distribution  is 

7 

6fb(z,Pj,t(z))  =  ef  — 

0 

to  first  order  in  the  radiation  fields,  where  the  solution  is  parametrized 
in  terms  of  the  axial  distance  from  the  start  of  the  interaction  region 

7 

(at  z  =  0),  and  t(z)  =  t  +J  dz'/v,(z')  is  the  sum  of  the  time  required 

U  -0  " 

for  an  electron  to  traverse  the  distance  and  the  entry  time  tg.  The 
transverse  components  of  the  trajectory,  therefore,  are  written  in  the 
form  px  =  pw  coskwz+Pxcosfi0t(2)-Py  sinSigtfz),  and  py  =  pw  sinkwz  + 

P^sinft  gt(z)  +  P  cos^0t(z). 

We  assume  plane  wave  solutions  of  the  form  exp(iajt)  and  choose 
to  work  with  scalar  and  vector  potentials  of  the  form 

5t(z,t)  =  y  6$ (z) exp  ( - iut)  +  c.c., 

and 

5^(3, t)  =  j  o^(z)exp( -iwt)  +  c.c., 


dz' 

TTzT 


(6|(z',t(z')) 


k ^ 


>5£(z',t(z')) 


W) 


(25) 


15 


where  it  is  evident  that  ' 5^  =  0.  After  transformation  to  the  basis 
i  i|y),  therefore,  the  perturbed  distribution  can  be  written 
as  '5fb(z,j3,T(z))  =  6f^(z ,^)exp(-iojt(z))  +  c.c.,  where 


7 


[-cp  (z')3,-6$(z")+iw(p_(z')6A+(z') 


+  exp (inQt (z ’ ))  (ioj-vz Cz')3z,)SA+(2 
+  exp(-iQ0t(z'))  (iu-vz(z')32 ,)cA_(z’) 


1  3 

P  9P 

/  3  4 

•  9  \ 

[w 

1 

v x 

y / 

,h  . 

■  i  3 

'  l3Px 

3P 

y 

J 


(24) 


/\  ✓v 

p+  =  px  ~  ipv,  5A+  si(5Ax  +  i6Av),and  t(z.z')  =J  dz"/vSz").  Using 

—  y  —  "  A  Z  /«. 

(24)  the  perturbed  current,  6^(z,t)  =  (6J+(z)^++  6J_(z)§_+6J,(z)e,)exp(-i(jjt) 

+  c.c.,  can  be  computed  as  follows 


roJ 


(5J 


*  ‘  I  /  dPxdPydp  y! jf  %-p>l 


P+ 


(25) 


Since  the  assumption  of  small  P  and  Py  is  central  to  the  analysis, 
we  shall  adopt  an  equilibrium  distribution  of  the  form 


Fb(px,pv,p)=  nbc(px)^(pv)Gb(p) 


(26) 


in  the  interest  of  computational  simplicity.  Here  G.(p)  is  an  arbitrar>r 
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function  of  the  magnitude  of  the  momentum  subject  only  to  the 
normalizationy*  dppG^(p)/pz  =  1.  Using  (26),  we  find  the  perturbed  currents 


5J 


i '  ^^”dp  V 


/  p+p.\  p + 

exp(+  i£2qt(z))  (2+ — —  P++  exp(+iQQT(z))  -r-  D_ 

\  p:  /  -  p;  + 


+  P+  (w~  +  hr)  D++p+(aF‘  -  HP")  D--P^Dz  Ip  | 

'  x  y/  \  x  v/ 

P  =P  =0 
x  v 


(27) 


6J 


z  ■  8 wf0  dP  t  (-IF  ' D» 


+  1  3P  _13P3  '  D 

x  v/ 


3p 


^(p),  (28) 


VPv"° 


where  pT  -  p  =  (p2-p„)?i,p+  =  p  exp(T  ik^) ,  t(z,z')  =  (z-z')/v 


U)l 


=  4-rre"nb/m, 


D+  =  -exp  (+ify  x  ( z) )  (5A+  ( z)  - <$A+  (0 )  exp  ( i  (u+f^ )  z/ vj 


;  i  dz"<SA+(z')exp(i(co+f20.lT(z,z' 


(29) 


and 


D.  =  -/  dz'exp(iajT(z,z')) 
u 


-9zMz')+  <SA+(z')+  —  6A_(z')/J  (30) 


P. 


P+ 


1 
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for  Px  *  P  =  0.  The  dispersion  equations  are  obtained  by  using  these 
perturbed  currents  in  the  wave  equations 

(>l  *  -pr)  4  «.• 

3,  mz)  -  sjz.  ^  ’ 

We  are  primarily  interested  in  the  low  gain,  tenuous  beam  limit 

in  which  <<  oj  and  the  wave  amplitudes  vary  little  over  the  length  of 

the  interaction  region.  This  is  the  regime  appropriate  to  the  experiments 

1  2 

conducted  at  Stanford  University.  ’  In  this  regime  collective  effects 
are  unimportant  and  the  space -charge  potential  can  be  neglected.  In 

addition,  the  coupling  between  the  electromagnetic  modes  described  by 

~  ->  7 

5A+  and  5A_  scales  with  P^/P_-  However,  the  assumption  of  low  gain  implies 
that  this  parameter  is  small,  and  we  can  focus  attention  on  one  or  the 

A 

other  of  these  modes  exclusively.  We  choose  to  consider  6A+,  which  in  the 
low  gain  regime  can  be  represented  as 


SA+(2)  =  6A+(0)exp(i^Z  dzTe+(z'))  ,  (32) 

where  |  link  (z)  |  «  |Rek+(z)|.  Under  the  assumption  that  k  (:)  =  k  +  5k+(;), 
where  k  is  independent  of  axial  position  and  |6k+(z)/k|  <<  1,  the  dispersion 
equation  becomes 


i"  im  win... i ii 

-  -»►. .  «i  wtk  ., 


k" 


2,  2 

u>  /c 


^  0, 


(35) 


and 


%  f°°  1  I  7  sin (.<*)/ v  - 

Im6k  (z)  %  -  - 7  /  dp  T  mw6“  - 7 - - U— 

*  -  2kc2 *  <  I2  K  “/v|l  k" 


sin(o)/v„-k-kw)z  8 


ii  R  K 

II  w 


+  JL 

YP, 


W 


0) 


v,  -n  ^ 

W  ||  0 


9P 


orfl 


w-(k+k nv.,11  k  v.  -P 

w  llV  w  0 


JL 


+  VifV(-vkv 


■] 


sin  (w/y,  -k-k  ) z 
w 


yp7 


ftttAI  „.o  -kv..  sin(“/vifVli  ‘k)z 


W  0 


co-kv, 


(34) 


o  “’ll 


II  0'  II 


Vp) 


1  ~>  7 

where  8“  =  p“/p“.  This  result  reduces  to  that  found  by  Sprangle  and 

W  w  z> 

4 

Smith  in  the  limit  of  a  vanishing  axial  magnetic  field. 

The  total  gain  over  an  interaction  region  of  length  L  is  defined 
as  G.  =  - f k  dz'Im<Sk+(z) .  Integrating  (34),  therefore,  we  find 


“b1  rZr  r^\^l  9  fn-  Sm  °w 

'  -7^XdpGb(P)  ~  3 pls„-T- 


2kc“  *'0 


n  6“  w 
P  w 


sin“0. 

e4 

w 


!)L  P  \  sin"0 


*z  ^  V  +  2vl|  “®e )  0W  ■ T  YPz  ^  2v||  e- 


l  „  6  PJ  / 

W  +  _E_  w  JL.M 


8"  PL  /  u)  +  Pn  \  sin"9, . 
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_E_ 

YP. 


8‘ 

1  +  ^  I1  + 

\y\\  “  0 


PnL\  sin“0, 

0  +  -- 

e  2  v;,  i 
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0: 


(55) 
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1 


where  9  =  (w/vM-k-k  JL/2  and  0„  5  (oj/v, -Q  /v,.-k)L/2.  If  the  beam 

w  II  *  8  II  0  || 

is  sufficiently  cold  that  a  distribution  of  the  form 

P, 

^(P)  =  -j  5(P-P0)  (36) 


can  be  used,  then  the  total  gain  becomes 
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(57) 


1  +  3°  (1  + 


KV\\o'% 


1  - 


fioL 


sin 
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‘V|j  0  egO 
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*>  7  u 
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In  Eq.  (57)  >  Yq  =  (1  -  P^/m-c-)  2,  Y,0  =  d  -  vj/c V',  3;0  =  P^/P^, 

PwO  E  %  pr0/(fi0'kwv||0r’  and  P||0  ^  V||0(E  p| |  0//y Om)  are  determined  by 
the  appropriate  solution  to  Pq  =  p"q  +  pj,g.  It  is  clear  from  (37), 
therefore,  that  wave  amplification  can  occur,  in  principle,  for  frequencies 


corresponding  to  both  the  usual  free  electron  laser  resonance 

(oj  2(k+kw)Vji  q)  and  the  Doppler-shifted  gyroresonance  (u  ^  ^  +  kvj|  q)  . 
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If  the  pump  period  is  short  compared  to  the  length  of  the 

system,  the  gain  for  frequencies  oj  'C  (k+k  )v,  is  dominated  bv  the 

w  ||  0 

first  term  in  (37) 


2  T  \ 

■y  w,  L  k 

Gi  2  — X — r 

16y  y  v, 

0  ZQ  10 


Yz0«w0°0  \  3  /sin6wC 

^0(-1+ew0'l’kwV||0  '  30wo  '  9w0 


The  extrema  occur  for  &wq  “■  +  1.3  at  which  3(sin0  q/S^q)  V30  q  ^  +  *54, 


which  correspond  to  frequencies 


^oW1 1  2'6/k“L); 


and  peak  gains 
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t!0tl*ew0,'\y||0 


where  the  plus  or  minus  sign  must  be  chosen  to  give  a  positive  gain  in  (40) . 
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These  expressions  for  the  gain  are  similar  to  those  found  previously  ’ 
in  the  limit  of  zero  axial  field.  The  principal  diferences  between  this 
result  and  those  found  previously  are  that  (1)  gain  can  occur  for  both 
frequencies  shown  in  (59),  and  (2)  a  singularity  is  introduced  for  kwv =“ 


P.0(l  +  3^q )  which  implies  that  the  effect  of  the  axial  guide  field  can  be  to 


substantially  enhance  the  small  signal  gain  over  the  case  in  which  no  guide 
field  is  present.  Note,  however,  that  the  singularity  is  a  product  of  our 
choice  of  a  monoenergetic  beam.  The  inclusion  of  an  energy  spread  can  be 
expected  to  broaden  the  resonance  and  remove  the  singularity. 

We  now  focus  on  the  effect  of  the  axial  magnetic  field  on  the  gain 
for  the  free  electron  laser  resonance.  If  A  denotes  the  ratio  of  the  gain 
in  the  presence  of  an  axial  guide  field  to  the  appropriate  expression  in  the 
limit  BQ  -*■  0,  then  it  follows  that 

2  eV  +  (U'l)  [e2  +  (y-D2Yo2  j 

A  =  -  U  — n - -r- — 7 - —  »  (41) 

[e  -  (u-1)  j[e  +  (y-1)  Y0  J 

where  u  =  k  vi^/H-.,  and  v  n  must  be  computed  in  a  self-consistent  manner  to 

W  ||U  U  jjU 

recover  the  constant -v^solutions  in  (5).  As  a  result,  A  can  be  expressed  as 
a  function  of  e,  S  ,  and  g^  (where  the  choice  of  constant -v,  solutions  im- 

VO 

plicitlv  selects  8  ) .  With  these  considerations  in  mind  we  plot  A  versus 
3q  for  fixed  values  of  total  electron  energy  and  wiggler  period  and 
amplitude  in  Fig.  9c.  IVe  shall  restrict  consideration  to  stable  uniform-v- 
traiectories  propagating  parallel  to  the  axial  guide  field.  We  choose 
3^  =  .97  and  S^/k^.c  =  .07  and  plot  the  variation  in  Vj  |q  with  8q  in  Fig.  9a. 
Evidently  there  is  a  stable  high  axial  velocity  trajectory  (orbit  I)  in  the 
1ow-3q  regime  (<  .72  for  the  chosen  parameters),  as  well  3S  a  trajectory 
which  varies  widely  in  magnitude  with  Bq  (orbit  II).  Since  the  frequency 
of  the  amplified  modes  scale  approximately  as  v| |q»  an  electron  beam 
characterized  by  trajectories  in  group  I  can  be  expected  to  excite  high 
frequency  waxes  for  all  accessible  values  of  3_.  In  contrast,  trajectories 


in  group  II  can  give  rise  to  a  wide  spectrum  of  oscillations  and,  for  the 
chosen  parameters,  and  high  frequency  waves  will  result  only  for  8q  >  1. 

This  is  clearly  shown  in  Fig.  9b  in  which  we  plot  the  resonant  frequency 
versus  8q  for  both  classes  of  orbits .  The  significant  feature  which  we 
are  concerned  with,  however,  is  the  effect  of  the  axial  guide  field  on 
the  small-signal  gain.  To  this  end,  we  plot  A  versus  8q  in  Fig.  9c  for 
the  stated  parameters.  We  observe,  first,  that  A  _>  1  throughout  the 
entire  range  of  8q  accessible  to  orbits  in  group  I  and,  hence,  the  gain 
is  enhanced  relative  to  the  zero  axial  field  (Sq  -*■  0)  limit.  In  addition, 
significant  enhancements  in  the  gun  are  possible  near  the  indicated  resonance 
at 


S 


0 


(1  + 


(42) 


which  is  the  orbital  stability  boundary  discussed  previously  in  the 
regime  in  which  Sq  <  v^/c.  The  resonance  condition  is  not  found  for 
orbits  in  group  II,  however,  since  sQ  >  v||o/c  f°r  aH  such  trajectories. 
Indeed,  A  vanishes  for  Sq  -  1.39  in  the  range  of  axial  fields  studied. 
Enhancement  in  the  gain  for  these  orbits  occurs  only  for  limited  ranges 
of  Sq.  A  significant  enhancement  is  found  in  the  1ow-Sq  regime  (i.e., 

•1  £  3q  <  1.34),  but  corresponds  to  comparatively  low  frequency  waves 
(see  Fig.  9b).  Enhancements  in  the  gain  for  high  frequency  waves  is 
found  only  for  a  limited  range  of  axial  field  strengths  (1.45  <  Sq  $  1.8), 
but  is  at  most  of  the  order  of  35%. 

The  enhancement  in  the  gain  for  axial  fields  such  that  Sq  is  of  the 
order  of  unitv  is  readily  explained  from  consideration  of  the  behavior  of 
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vw>  It  is  clear  that  for  axial  fields  in  which  fig  ^  kwvj|Q  (i.e.,  8g  ^  1) 
large  enhancements  in  the  wiggle  velocity  vw  occur.  This  is  the  source  in 
the  enhancements  in  both  the  gain  and  emissivity.  However,  as  Bn  continues 
to  increase  vw  -►  eVy  g  and  the  wiggle  velocity  decreases  relative  to  that 
found  in  the  zero  axial  field  limit. 

In  the  absence  of  an  axial  guide  field  the  small  signal  gain  scales 
as  the  square  of  the  wiggler  amplitude  (i.e.,  q^k^c  )  which  is  also  the 
square  of  the  wiggle  velocity.  However,  for  finite  values  of  the  axial 
guide  field  this  variation  is  more  complex  and  the  gain  does  not  increase 
monotonicallv  with  either  Q^/k^.c  or  the  wiggle  velocity.  This  situation 
is  illustrated  in  Fig.  10,  in  which  we  plot  the  peak  gain  (solid  line)  and 
frequency  corresponding  to  peak  gain  (dashed  line)  versus  q^/k^.c  for 
8 g  =  1.13  and  v/c  =  .96.  For  small  values  of  this  parameter  the  gain  does 
indeed  scale  as  the  square  of  the  wiggler  amplitude,  but  a  peak  occurs  at 
qykwc  -012  after  which  the  gain  rapidly  drops  to  zero.  In  addition  the 
frequency  monotonically  decreases  throughout.  Note  that  while  o/l^c  does 
not  equal  the  wiggle  velocity  in  the  presence  of  an  axial  field,  the 
decreasing  frequency  implies  a  decreasing  axial  velocity  which,  in  turn, 
leads  to  an  increasing  wiggle  velocity.  This  is  illustrated  more  clearly 
in  Fig.  11  in  which  we  plot  the  variation  in  the  axial  velocity  versus 
q;v/kwc  for  v/c  =  .96  and  Sg  =  .75  and  1.13.  For  8g  =  1.13  all  orbits  are 
of  group  II  class  and  the  axial  velocity  is  a  decreasing  function  of 
“w/^w c  over  t^ie  ran8e  studied.  Since  this  means  that  the  wiggle  velocity 
increases  with  the  wiggler  amplitude  over  this  range,  it  must  be  concluded 
that  the  small  signal  gain  cannot  necessarily  be  increased  by  a  simple 
enhancement  in  the  wiggle  velocity. 
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For  sufficiently  smaller  values  of  Bg,  group  I  orbits  appear  (see 
Fig.  11  for  Sq  =  .75)  over  a  limited  range  of  wiggler  field  amplitudes. 

The  peak  gain  (as  well  as  the  corresponding  frequency)  for  this  case  is 
shown  in  Fig.  12,  and  describes  the  frequency  and  gain  associated  with  the 
group  I  orbits.  We  have  restricted  consideration  to  group  I  orbits,  because 
while  group  II  orbits  also  lead  to  growth,  the  gain  is  weaker  and  the 
frequency  is  relatively  low  (w  £  8  k^c) .  For  these  orbits ,  the  gain  is 
seen  to  increase  with  increasing  wiggle  velocity,  and  the  singularity 
mentioned  previously  is  also  found  at  the  transition  to  orbital  instability. 
It  should  be  remarked,  however,  that  the  small  signal  approximation  breaks 
down  in  the  vicinity  of  the  singularity. 

The  gain  for  gyroresonant  emission  is  given  by  («  fig  +  kvj.g) 


for  a  long  interaction  region  (kl,  >>  1) ,  and  growth  is  possible  if 

Kv\\o  *  Qo 


8w0 


s  fro  -  Vhq) 

*  Ck  ♦  Is/llo 


(44) 


Peak  amplification  occurs  at  the  maximum  of  (sine^/e^) "  1  1,  which 
corresponds  to  a  frequency  u  =  ftp  +  kvng  and  peak  gain 


i  1 
(Ur  L 

^GlWx  -  T 


2y0c 


v||0/c^20fl0  +  V 

°o""  kwV||  0 


(45) 
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The  variation  of  the  peak  gain  and  associated  frequency  for 
gvroresonant  emission  with  axial  field  strength  are  shown  in  Fig.  13  for 
parameters  consistent  with  the  calculation  of  the  frequency  and  peak  gain 
at  the  free  electron  laser  resonance  presented  in  Fig.  9.  Gyroresonant 
emission  is  found  only  for  group  II  orbits  and,  as  expected,  the  frequency 
increases  with  increasing  8q.  In  contrast,  the  gain  is  found  to  decrease 
with  increasing  axial  field  strengths.  The  large  increase  in  the  gain  as 
8q  approaches  zero  results  from  the  fact  that  Vjj  q  also  vanishes  in  this 
limit.  However,  the  small  signal  gain  approximation  breaks  down  in  this 
limit,  and  nonlinear  effects  can  be  expected  to  produce  enhancements  lower 
than  expected  on  the  basis  of  the  linearized  theory. 
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V. 


SUMMARY  AND  DISCUSSION 


In  this  work,  we  considered  both  the  spontaneous  and  induced 
emission  from  a  free  electron  laser  in  which  a  uniform  axial  guide  field 
is  present.  The  existence  of  the  axial  field  is  found  to  introduce  an 
additional  source  of  coupling  between  the  axial  and  transverse  electron 
motion  which,  for  sufficiently  strong  field  levels  (’I),  can  lead  to 
orbital  instability  for  motion  parallel  to  Since  the  existence  of 
orbital  instability  results  in  a  rapid  dispersal  of  an  initially  bunched 
electron  beam  with  disastrous  implications  for  the  expected  radiation 
levels,  we  focus  attention  on  the  parameter  regimes  corresponding  to 
orbital  stability  and  orbits  characterized  by  nearly  uniform  axial  velo¬ 
cities.  The  predominant  transverse  motion  of  these  orbits  is  to  track 
the  wiggler  field;  however,  a  relatively  small  component  of  Larmor  ro¬ 
tation  associated  with  the  guide  field  may  also  be  present.  Because  of 
this ,  emission  can  be  expected  for  both  the  free  electron  laser  and 
gyrotron  resonance  conditions. 

The  spontaneous  emission  describes  the  ambient  level  of  noise 
within  the  device  prior  to  the  coherent  amplification  process.  .As  shown 
in  Eq.  (21) ,  the  spontaneous  emissivity  for  waves  propagating  axially 
contains  components  due  to  the  transverse  motion  resulting  from  the 
action  of  the  wriggler  and  axial  guide  fields.  In  each  case  the  level 
of  emission  is  proportional  to  the  square  of  the  magnitude  of  the  trans¬ 
verse  velocities.  This  velocity  corresponds  to  the  canonical  momentum 
in  the  transverse  direction  for  the  gyrotron  mode  (u  =  fig  +  kv||  ) .  and 
to  the  wilier  velocity  v  for  the  free  electron  laser  mode 

W 


(u  =“  (k+kw)v  ) .  It  should  be  recognized  that  the  presence  of  a  finite 
leads  to  an  enhancement  in  v  and,  consequently,  in  the  free  electron 
laser  emissivitv  bv  a  factor  of  (ftn/k  v  -  1)  .  As  discussed  in  Sec.  IV 

4  '  U  \v  j; 

the  limit  in  which  1  =  k  v  is  not  achievable. 

o  w  „ 

Enhancement  in  the  linearized  (or  small  signal)  gain  also 
expected  to  result  when  an  axial  guide  field  is  present.  For  the  ten¬ 
uous  beam,  low  gain  regime  the  enhancement  of  the  gain  is  expected  to  be 
greatest  for  parameters  in  the  vicinity  of  the  transition  to  orbital 
instability  (denoted  by  the  dashed  line  in  Fig.  9)  which  corresponds  to  a 
beam  characterized  by  orbital  stability.  In  order  for  the  analysis  to  be 
valid  for  orbitally  stable  beams,  however,  stringent  requirements  on  beam 
quality  must  be  satisfied  because  a  small  energy  spread  can  lead  to  a 
rapid  breakdown  in  the  coherence  of  the  beam.  In  addition,  the  small - 
signal  gain  approximation  must  break  down  in  the  vicinity  of  the  resonance 
(i.e.,  the  transition  to  orbital  instability) .  Finally,  significant 
enhancements  in  the  gain  are  expected  to  occur  for  low  wavelength 
oscillations  only  in  the  low-gg  regime. 
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Fig.  5  -  Graph  of  the  bounds  on  the  axial  velocity  versus  g 

v/c  =  .97,  and  6q  =  .1. 
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Fig.  9  -  Graphs  of  axial  velocity,  and  resonant  frequencies  and 
enhancements  in  the  gain  for  stable  trajectories  versus  6q  for 


v/c  =  .97  and  a 7k  c  =  .07. 


.01  .02 
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Fig.  10  -  Graphs  of  the  maximum  gain  (solid  line)  and  the  corresponding 
frequency  (dashed  line)  versus  wiggler  strength  for  =  1.13 
and  v/c  *  .96. 
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